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SYMBOLS 


D/,Sw 

aircraft  drag  coefficient 


local  wing  chord 
maximum  store  diameter 
store  length 
Mach  number 

free  stream  dynamic  pressure,  pof 

2 

8 tore  msx  cross  sectional  area  nd  /4 
pylon  planform  area 
store  tail  wetted  area 
aircraft  wing  area 

lateral  distance  from  store  center  line  to  aircraft 
center  line 

minimum  depth 

maximum  depth,  (see  Figure  3) 

equivalent  body  length  (see  Figure  4) 

wing  leading  edge  Bweep  angle 

distance  from  store  nose  to  wing  leading  edge 

equivalent  body  diameter  (see  Figure  4) 
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SUMMARY 


The  drag  due  to  the  external  carriage  of  weapons  on  high  speed 
attack  aircraft  results  in  large  reduction  in  speed  and  range.  Up 
to  the  present,  the  only  means  of  evaluating  this  carriage  drag  has 
been  either  full  scale  flight  test,  or  wind  tunnel  testing.  Both 
methods  are  costly  and  time  consuming.  Accordingly,  an  analytical 
method  is  being  sought  which  will  provide  useful  estimates  of  the 
drag  increment  caused  by  external  weapon  installation  without  the 
need  for  such  tests.  The  method,  based  on  regression  techniques, 
utilizes  the  large  amount  of  wind  tunnel  data  obtained  from  previous 
tests  at  NSRDC  and  other  facilities.  The  prediction  equation  is  de¬ 
termined  as  a  function  of  the  significant  geometric  parameters  of 
the  aircraft-store  configuration,  A  highly  encouraging  degree  of 
success  has  been  achieved  in  application  of  the  method  to  drag  of 
stores  mounted  singly  on  pylons  at  Mach  numbers  of  0.8  and  0.9. 

The  prediction  equation  has  shown  excellent  agreement  with  drag 
data  subsequently  measured  on  arbitrary  aircraft-store  combinations. 
Effort  to  extend  the  scope  of  application  is  continuing. 

INTRODUCTION 

The  drag  caused  by  externally  mounted  weapons  or  Captive  Flight 
drag  has  been  a  problem  area  in  Aerodynamics  for  nearly  as  many  years 
as  man  has  been  flying  aircraft.  During  World  War  II,  while  there  were 
performance  losses  due  to  external  carriage  of  weapons.  fho«e  lessee 
were  not  sufficiently  large  to  cause  alarm.  As  speeds  of  aircraft  in¬ 
creased,  the  problem  became  more  critical.  Whereas  a  World  War  II 
aircraft  might  auffer  a  reduction  in  speed  of  20  miles  an  hour  due  to 
external  fuel  tanks,  todays  attack  aircraft  may  lose  100  to  200  MPH, 

An  example  of  this  type  of  loss  is  shown  in  Figure  1.  It  might  well 
be  pointed  out,  that  the  situation  depicted  is  not  an  extreme  case  and 
that  there  are  many  worse. 


The  initial  investigations  into  the  matter  of  captive  flight  drag 
at  N3RDC  were  of  limited  scope.  These  investigations  were^  in  general, 
merely  determining  experimentally  the  drag  of  a  particular  weapon 
carried  on  a  particular  aircraft.  While  these  investigations  pre¬ 
sented  useful  data,  they  presented  very  little  insight  into  the 
problems  involved.  It  was  not  until  1960  that  a  systematic  investi¬ 
gation  was  made  into  the  overall  problem  of  captive  flight  drag. 

Initial  tests  involved  research  models  known  as  MIDS  I,  II,  and  III. 

The  variables  studied  included  wing  sweep,  store  position,  wing  leading 
edge  droop,  nose  shape,  and  others.  At  the  same  time  that  these 
studies  were  being  conducted,  additional  testing  involving  current 
day  Navy  aircraft  was  initiated.  Aircraft  such  as  the  A-4,  A-7A, 
and  F-8  were  tested  with  many  of  the  possible  weapon  loadings  and  the 
results  published.  Limited  wind  tunnel  model  investigations  were 
also  carried  out  on  several  Navy  aircraft  in  order  to  reduce  the 
drag  of  stores,  and  set  up  guidelines  for  future  attack  aircraft. 

By  1966  sufficient  experimental  data  had  been  accumulated 

from  NSRDC,  NASA,  and  other  sources  to  warrant  assembly  into  a 
convenient  and  usable  form.  Accordingly,  data  for  aircraft  with  and 
without  external  weapons  were  recorded  on  punched  cards  and  a  computer 
program  written.  This  program  curve-fitted  and  solved  for  the 
incremental  drag  values,  and  automatically  plotted  the  results  to  a 
suitable  scale. 

Very  little  attempt  at  correlation  of  these  available  data  with 
significant  parameters  of  the  aircraft-store  combination  had  been 
made  In  the  past.  It  was  decided,  that  by  using  the  high  speed 
digital  computer  and  statistical  techniques,  correlation  might 
successfully  be  attempted. 
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BACKGROUND  FOR  METHOD  OF  SOLUTION 
It  is  well  established  from  both  experimental  and  theoretical 
aerodynamics  that,  for  practical  engineering  purposes 

CD  *  CD^®  *  Re*  Body  Geometry) 

Within  the  ranges  of  variables  of  greatest  immediate  interest,  differ¬ 
ences  in  incremental  drag  due  to  Reynolds  number  are  usually  of  secondary 
importance,  and  it  was  decided  to  neglect  these  differences. 

The  key  to  successful  correlation  is  to  describe  the  body  geometry 
in  terms  of  a  manageable  number  of  the  most  important  nondimens ional 
geometric  parameters.  These  might  include  nondimens ional  combinations 
of  length,  cross  sectional  area,  diameter,  etc.  The  drag  coefficient  at 
a  given  Mach  number  might  then  be  written  as 

CD  =G(V  V  *  •  •  *  Xn) 

where  x^,  .  .  .  ,  x^  are  geometric  parameters. 

If  is  thought  to  be  a  linear  expansion  of  the  parameters,  then 

Ch  =  a  +  a.  x.  a.  x„  +  ,  .  ,  a  x 
Doll22  nn 

Additionally  it  might  be  necessary  to  use  higher  order  powers  of  the 

parameters  and/or  cross  product  terms,  so  that 

C_  =  a  +  a,  x-  +,...+  a  x  +  b.  x2,  +...,+  b  x  2 
Doll*  n  n  11  n  n 

"b  X-  Xn  *b  •  •  •  *b  C  1  X  1  X  +  •  •  • 

12  1  2  n-1,  n  n-1  n 

For  the  kth  body  of  a  set  of  many  bodies,  the  equation  becomes 

CD,  S  ao  +  ai  Xik  +  bi  Xik2  +  Cij  Xik  Xjk  +  *  '  • 
k 

If  sufficient  experimental  data  are  available  the  unknown  coefficients 

a  »  a.,  b . ,  c . .  may  be  found  by  statistical  methods, 
o’  i*  i’  ij  J 

MULTIPLE  LINEAR  REGRESSION 

The  statistical  method  chosen  to  solve  for  the  unknown  coefficients 
is  multiple  linear  regression.  Multiple  linear  regression  involves  the 
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fitting  of  a  surface  in  a  "K"  dimensional  space.  This  surface  is  the 
best  possible  fit  of  an  assumed  value  such  that,  its  deviation  from  the 
true  value  is  minimized.  The  form  of  the  regression  equation  is 


a  +  a,  B 
o  1  w 


lk  +  a2  p2k  +  *  •  +  an  Pnk 


where  o is  the  estimated  value  of  and  the  j3^>  i  *  1,  n  are  the 
known  parameters. 

Measuring  the  parameters  from  their  means,  the  equation  may  be 
rewritten  as 


V'  -  V'  hk  +  V  ®2k"  +  •  •  •  +  a„"  *nk" 


where 


°ic  "  \  “ 


and 


*ik"  -  plk'  •  h>  -  ai 


and 


ao  “  al  ^  +  • 


+  an  K  - 


it  can  be  shown  (see  the  appendix)  that,  from  the  above,  a  system  of 
equations  may  be  generated  with  respect  to  the  least  squares  criterion, 
and  the  unknown  a^*s  may  be  found.  The  regression  equation  has  the  same 
form  as  the  equation  for  if: 

V  -  \ 

ao  -  V  ai  *  V  bi  or  cij 


and 


ik 


ik*  xik 


or  xik  xjk 


In  addition  to  the  solution  of  the  unknown  a^'s,  two  additional  coeffi¬ 
cients  may  be  obtained  from  the  regression  equation.  The  first  is  called 
the  multi-correlation  coefficient.  This  coefficient  gives  a~  indication 
of  the  goodness  of  fit  of  the  regression  equation.  The  second  coefficient 
is  called  the  partial  correlation  coefficient.  This  indicates  the  impor¬ 
tance  of  a  single  parameter  in  the  overall  equation.  The  mathematical 
definition  of  these  coefficients  is  given  in  the  appendix. 
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PARAMETERS 


The  parameters  which  are  used  tu  describe  the  geometry  of  the 
wing  pylon-store  configuration  include  store  position  relative  to  the 
wing,  reference  areas  of  store  and  pylon,  store  fineness  ratio  and  wing 
leading  edge  sweep  angle.  In  addition,  at  certain  Mach  numbers,  area 
rule  effects  become  important  as  well  as  drag  interference  between  store 
upper  surface  and  wing  lower  surface. 

With  the  inclusion  of  an  area  rule  parameter  and  interference 
parameters  the  final  correlation  equation  for  application  to  data  at  a 
single  Mach  number  is: 

ACD  “  ao  +  al  (VO  +  a2  (1l/d)  +  a3  (?/c)  +  a4  (£/d) 

+  a,  (s/S  )  +  a,  (1/d;  +  a,  (s„/S  'j  +  a„  (s  /S  N 
5  V  p  w.'  6  v  '  7  \  t  m  '  8Vmwy 

+  ag  (cos  X)  +  a1Q  ((t/h)cos  X)  +  an  (V^x)3 

+  a12  (Sm/Sw)2  +  a13  ^cos  +  a14  ((T/H)cos  +  a15  (¥d)3 

+  °16  (^U)  Wd)  +  a17  (^/TUx)  (5/c> 

The  various  parameters  are  shown  in  Figure  2. 

The  interference  parameter  is  shown  in  Figure  3.  The 

ratio  approximates  the  familiar  A/A*  of  Laval  nozzle  theory.  When  the 
ratio  becomes  small,  the  chance  of  sonic  flow  becomes  large  and  therefore, 
large  interference  may  occur. 

Limited  area  rule  has  been  taken  into  account  by  the  use  of  the 
parameter  ((t/h)cos  X)  >  and  is  shown  in  Figure  4.  The  ratio  (t/h)  *  is 
the  fineness  ratio  of  an  equivalent  body  encompassing  part  of  the  wing, 
pylon,  and  body  as  shown  in  Figure  4. 

Using  the  correlation  equation,  good  agreement  between  predicted 
and  experimental  data  has  been  obtained.  This  is  shown  in  Figure  5  for 
Mach  numbers  of  0.8  and  0.9.  Figure  5  presents  data  for  forty-nine  air¬ 
craft  store  configuiations  consisting  of  eleven  aircraft  and  fourteen 
stores.  The  data  are  presented  in  the  form  of  and  are  based  on  an 
accuracy  of  ±0.0010. 
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ASSESSMENT  OF  PARAMETERS 

The  importance  of:  various  parameters  on  incremental  drag  may  best 
be  assessed  by  use  of  the  aforementioned  partial  correlation  coefficient. 
The  effect  of  a  change  in  incremental  drag  caused  by  a  change  in  a  param¬ 
eter  is  determined  by  the  sign  of  its  coefficient  in  the  regression  equa¬ 
tion.  Utilizing  for  each  parameter,  the  partial  correlation  coefficient 
and  the  regression  coefficient  associated  with  the  parameter,  the  following 
were  observed: 

1.  The  interference  parameters  C Hnax )  aricJ  (^^lax)2  had  little 
effect  at  M  s  0.8.  At  M  *  0.9  interference  became  important  and  incre¬ 
mental  drag  decreased  with  an  increase  in  (T)/’f\nax)  •  *n  general  (j\^max) 
is  less  than  1  and  the  ratio  should  be  as  close  to  1  as  possible  for  low 
incremental  drag  at:  high  subsonic  Mach  numbers. 

2.  The  depth  beneath  the  wing  (T]/d)  and  (Ti/d)2  were  found  to  be  of 
equal  importance  at  both  Mach  numbers  evaluated.  An  increase  in  (7]/d) 
produced  a  corresponding  increase  in  incremental  drag. 

3.  The  longitudinal  position  parameter  (§/c)  was  important  at  both 
Mach  numbers.  The  incremental  drag  increased  with  an  increase  in  (£/c). 

4.  Variation  of  incremental  drag  with  lateral  position  (£/d)  was 
small.  An  increase  in  (£/d)  produced  a  smill  decrease  in  incremental  drag. 

5.  Pylon  size  in  relation  to  wing  area  was  important  at 

M  *  0.8.  Where  an  increase  in  (s  /S  )  produced  a  larger  incremental  dLwg. 
At  M  *  0.9  the  importance  of  the  parameter  was  not  as  noticeable  and  the 
effect  on  drag  was  reversed. 

6.  The  size  of  the  stores  tail  surfaces  in  relation  to  cross  sec¬ 
tional  area  (s^/S^was  important  at  both  Mach  numbers.  The  effect  on 
drar>  was  constant  at  both  M  *  0.8  and  M  *  0.9.  As  the  ratio  increased, 
the  incremental  drag  increased. 

7.  The  effect  on  drag  caused  by  a  change  in  fineness  ratio  (1/d) 
was  important  at  both  M  *  0.8  and  M  *  0.9.  As  expected  an  increase  in 
fineness  ratio  produced  a  decrease  in  incremental  drag. 

8.  The  size  of  a  store  carried  on  a  particular  aircraft  (sm/sw)> 
and  ( S  /S  )s  was  of  relatively  minor  importance  at  M  «  0.8.  At  M  *  0.9, 
however,  the  effect  was  more  pronounced.  An  increase  in  ^ S^/S^ produced 
a  corresponding  increase  in  incremental  drag. 
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9.  Wing  leading  edge  sweep  angle  (cos  X)  and  (cos  X)2  was  of 
relative  importance  at  both  Mach  numbers.  At  M  *  0,8,  an  increase  in 
cos  X  ,  produced  a  decrease  in  incremental  drag.  At  M  *  0,9  the  effect 
on  drag  caused  by  an  increase  in  (cos  X)  was  reversed. 

10.  At  M  «  0.8  the  limited  area  rule  parameters  ((t/h)cos  X)  and 
((t/k)cos  X)2  have  little  effect.  At  a  Mach  number  of  0.9  the  parameter 
was  of  large  importance.  An  increase  in  the  parameter  produced  an  in¬ 
crease  in  drag. 

11.  The  cross  product  term  ( VUnav )(7]/d)  was  of  extreme  importance 
at  both  Mach  numbers.  The  effect  on  drag  was  opposite  to  that  produced 
by  (T|/ d )  and  (Tj/ d)*5.  An  increase  in  the  cross  product  parameter  produced 
a  decrease  in  incremental  drag  at  both  Mach  numbers. 

12.  The  other  cross  product  parameter  (<n/^\nax)  (£/c  was,  likewise, 
of  large  importance  at  both  Mach  numbers.  An  increase  in  the  parameter 
produced  a  decrease  in  incremental  drag. 

FUTURE  PLANS 

At  present  the  prediction  works  well  for  singly  mounted  bodies  of 
revolution  in  the  transonic  speed  range.  Attack  aircraft,  however,  are 
not  limited  to  a  single  carriage  per  wing,  but  often  have  provisions  for 
multiple  carriage  of  stores.  Many  of  the  stores  carried  are  not  plain 
bodies  of  revolution  but  have  complex  shapes,  such  as  multiple  bomb  racks 
and  winged  missiles.  Plans  are  being  made  to  incorporate  the  effects 
of  these  complex  configurations  in  the  prediction  equation.  Multiple 
bomb  racks  may  require  the  use  of  a  fictitious  body  having  the  same  drag 
characteristics  as  the  rack,  as  well  as  additional  parameters  for  proper 
drag  determination.  The  problem  of  side-by-side  store  mounting  will 
undoubtedly  require  additional  terras  in  the  equation  in  order  to  account 
for  store-to-store  interference. 

Optimization  of  store  position  for  minimum  drag  will  also  be 
attempted.  Wind  tunnel  tests  at  Naval  Ship  Research  and  Development 
Center  indicated  that  significant  gains  in  performance  can  be  attained, 
by  proper  store  positioning.  Use  of  the  correlation  equation  for 
purposes  of  optimization  is  very  promising,  and  does  not  require  the  time 
and  expense  associated  with  wind  tunnel  testing. 


As  more  data  become  available  the  equation  will  be  further  refined 
and  its  speed  range  extended  to  low  supersonic  speeds. 


Aerodynamics  Laboratory 

Naval  Ship  Research  and  Development  Center 

Washington,  D.  C.  20007 

July  1968 
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APPENDIX 

MULTIPLE  LINEAR  REGRESSION 

Given  n  sets  of  data  in  terms  of  k  parameters,  suppose  that  a 
dependent  variable  or  yield  parameter,  X^,  may  be  approximated  by  a 
linear  combination  of  the  other  k  -  1  parameters.  Then  a  general  equa¬ 
tion  would  be  of  the  form 

V  “  Bo  +  52  *2  +  •  •  *  Bk  \ 

where  X^  7  is  the  estimated  value  of  X^,  and  Xj,  j  *  2  .  .  .  k  are  the 
k  -  1  parameters. 

Measuring  the  parameters  from  their  means,  the  equation  may  be 
rewritten 

V  "  bo  +  b2  x2  +  ‘  •  •  bk  *k 

where 

*  X^  —  j  i  *  1,  •  .  .k 

bo  “  B2  *2  +  •  •  •  Bk*k  +  Bo  '  X1 
bj,  -  Bi,  i  -  2,  .  .  .  k 

For  X^ 7  to  be  the  best  linear  estimate  of  X^,  then  with  respect 
to  the  least  squares  criterion,  -  X^  7  'j  3  must  be  minimized.  This 

is  equivalent  to  minimizing  ^  “  x^2  ;  i.e., 

Ik  -  b0  -  b2  x2  -  •  •  •  bk  \T 

A  necessary  and  sufficient  condition  for  minimization  of  this  is 

I  I  £-  K  -  bo  *  b2  *2  -  •  •  ■  bk  *k)a  -  0 

j-o  J 

J# 
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writing  out  the  k  equations  gives 


b-“*fc*I*i*  -  -  -  sF-k  ■  I«, 

b.  5\  “.I' >*«»♦•  • 

Hole  -n*  £  “  £  (xt  -  i1)  -  0.  Ihut  etch  tera  tn  the  (ir,t  equation 

vanishes ,  and  bQ  -  0.  The  system  becomes 


\h‘*- 

.  .  .tbk^X) 

b2  Jx2  xk+  •  •  .  +  bk  £*2a  -Yxlxk 


These  are  the  normal  equations 
The  coefficients  of  the  b 
statistical  quantities: 


,  solvable  for  i  ri  2,  .  .  yit 
’s  can  be  expressed  in  the  following 


r 


ij 


sample  standard  deviation  of 

Ixi  xj 


-  sample  correlation  coefficient  between 

Xi  and  xj 

Then  r-* 

Zxi  xj  ■ n  rij  si  sj 
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Cancelling  common  factors  after  these  substitutions  have  been 
made,  the  normal  equations  become 


bZ  r22  *2  +  b3  r23  *3  + 
b2  r32  *2  +  b3  r33  *3  + 


*  ’  bk  r2k  “k 

•  *  bk  r3k  “k 


"  r21  *1 
"  r31  81 


b2  rk2  “2  +  b3  rk3  *3  + 


•  •  bk  rkk  8k 


rkl  81 


Assuming  that  the  determinant  of  the  coefficients  is  not  zero. 


81  82  •"  8i-l  “i+l  *•’  *k 


*2  *3  *"  8k 


r22 

r23  *  * 

*  r21 

r2k 

r32 

r33  ‘  * 
• 

• 

.  r31 

r3k 

rk2 

e 

rk3  ‘  ' 

’  rkl 

rkk 

r22 

r23  ’  ' 

'  "  r2i 

r2k 

r32 

r33  *  ' 
• 

• 

.  .  r31 

r3k 

rk2 

• 

r't3  *  ' 

'  *  rki 

rkk 

Note  that  the  determinants  in  the  numerator  and  denominator 
differ  only  in  the  elements  occurring  in  the  i  -  1st  column.  Con¬ 
sider  the  following  determinant: 


rll 

r12 

r13  * 

•  •  rik 

R  -  r21 

r22 

• 

ft 

r23  * 

•  *  r2k 

rkl 

• 

rk2 

rk3  * 

’  ’  rkk 
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The  determinant  in  the  denominator  of  is  the  minor  of  r^  in  R, 
denoted  R^.  Also,  the  numerator  determinant  of  can  be  converted  to 

the  minor  of  r^  in  R  by  shifting  the  column  headed  by  r^^  to  the  first 
column  position.  Using  the  relation  between  minors  and  cofactors,  the 
numerator  determinant  becomes  -  R^.  Thus 


b 


i 


The  usefulness  of  the  first  equation  for  estimation,  i.e., 


Xj '  -  Bo  +  B2  X2  +  .  .  .  +  ^  Xk 

is  indicated  by  its  multiple  correlation  coefficient  defined  in  Refer¬ 
ence  1,  as 

r1.23...k  " 

If  r^  k  is  dose  to  1,  this  indicates  that  the  data  points  lie 

near  the  derived  "plane,"  while  if  r^  23  ^  is  dose  to  zero  then 

either  the  relationship  is  weakly  linear,  or  curvilinear. 

Similarly,  the  influence  of  a  particular  parameter  is  reflected  in 
its  partial  correlation  coefficient: 


rlj.23...k 


-R 


11 


/hi  Rjj 
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Figure  8  -  Definition  of  Interference  Parameter  (rj/ijmmx) 
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Figure  5  -  Results  of  Prediction  Equation  for  49  Configurations 
at  Hach  Numbers  of  0*8  and  0.9 
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